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Various versions of the definition of nonclassical symmetries existing in the lit- 
' erature are analyzed. Comparing properties of Lie and nonclassical symmetries 

leads to the conclusion that in fact a nonclassical symmetry is not a symmetry 
r*j | in the usual sense. Hence the term "reduction operator" is suggested instead of 

Oh the name "operator of nonclassical symmetries" . It is shown that in contrast 

to the case of single partial differential equations a satisfactory definition of 
nonclassical symmetries for systems of such equations has not been proposed 
up to now. Moreover, the cardinality of essential nonclassical symmetries is 
discussed, taking into account equivalence relations on the entire set of non- 
classical symmetries. 
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■ 1 Introduction 

■ The "nonclassical" method of finding similarity solutions was introduced by Blu- 
man and Cole in 1969 [3]. In fact, the method was first appeared in [2] in terms of 
"nonclassical group" but the terminology was changed in [3]. Over the years the 
"nonclassical" method began to be associated with the term nonclassical symme- 
try [13] (also called Q-conditional [8] or, simply, conditional symmetry [6j [TO] ) . 
In the past two decades, the theoretical background of nonclassical symmetry was 
intensively investigated and nonclassical symmetry techniques were effectively ap- 
plied to finding exact solutions of many partial differential equations arising in 
physics, biology, financial mathematics, etc. See, e.g., the review on investigations 
of nonclassical symmetries in |18j . 

Here we mention only works which are directly connected with the subject of 
our paper. In the pioneering paper [3] the "nonclassical" method was described 
by means of the example of the (1 + l)-dimensional linear heat equation. It was 
emphasized that any solution of the corresponding (nonlinear) determining equa- 
tions gives the coefficients of an operator such that an ansatz based on it reduces 
the heat equation to an ordinary differential equation. A veritable surge of inter- 
est in nonclassical symmetry was triggered by the papers [161 \TT\ |9]. In |16| the 
"nonclassical" method was considered in the course of a comprehensive analysis 
of a wide range of methods for constructing exact solutions. The concept of weak 
symmetry of a system of partial differential equations, generalizing the "nonclas- 
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sical" method, was introduced in [T7], where also the reduction procedure was 
discussed. Moreover, fundamental identities |171 eq. (23)] crucially important for 
the theory of nonclassical symmetries were derived (see Myth [3] below). The first 
version of the conditional invariance criterion explicitly taking into account dif- 
ferential consequences was proposed in [S] . Generalizing results of [SJ [7] and other 
previous papers, in [6] Fushchych introduced the notion of general conditional 
invariance. From the collection of papers containing [6] it becomes apparent that 
around this time a number of authors began to regularly use the terms "condi- 
tional invariance" and "Q-conditional invariance" in connection with the method 
of Bluman and Cole. The direct (ansatz) method closely related to this method 
was explicitly formulated in [4]. To the best of our knowledge, the name "nonclas- 
sical symmetry" was first used in |13j . Before this, there was no special name for 
operators calculated in this approach and the existing terminology on the subject 
emphasized characteristics of the method or invariance. The involution condition 
for families of operators was first considered in the formulation of the conditional 
invariance criterion in jlO|, I27| . The relations between nonclassical symmetries, re- 
duction and formal compatibility of the combined system consisting of the initial 
equation and the invariant surface equation were discovered in [23] and were also 
studied in [15]. The problem of the algorithmization of calculating nonclassical 
symmetries was posed in [5]. Furthermore, the equivalence of the non-classical 
(conditional symmetry) and direct (ansatz) approaches to the reduction of partial 
differential equations was established in general form in |28j . making use of the 
precise definition of reduction of differential equations. 

In spite of the long history of nonclassical symmetry and the encouraging 
results in its applications, a number of basic problems of this theory are still 
open. Moreover, there exists a variety of non-rigorous definitions of related key 
notions and heuristic results on fundamental properties of nonclassical symmetry 
in the literature, which are used up to now and form what we would like to call 
the "mythology" of nonclassical symmetry. These definitions and results require 
particular care and presuppose the tacit assumption of a number of conventions in 
order to correctly apply them. Otherwise, certain contradictions and inaccurate 
statements may be obtained. Note that mythology interpreted in the above sense 
is an unavoidable and necessary step in the development of any subject. 

Basic myths on nonclassical symmetries presented in the literature are dis- 
cussed in this paper. We try to answer, in particular, the following questions. 

• Is a nonclassical symmetry a Lie symmetry of the united system of the initial 
equation and the corresponding invariant surface condition? Can a nonclas- 
sical symmetry be viewed as a conditional symmetry of the initial equation 
if the corresponding invariant surface condition is taken as the additional 
constraint? Is nonclassical symmetry a kind of symmetry in general? Does 
there exist a more appropriate name for this notion? 

• What is a rigorous definition of nonclassical symmetry for systems of differ- 
ential equations? Can such a definition be formulated as a straightforward 
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extension of the definition of nonclassical symmetry for single partial differ- 
ential equations? 

• Is the number of nonclassical symmetries essentially greater than the number 
of classical symmetries? 

2 Definition of nonclassical symmetry 

Following O [TOl [2B] , in this section we briefly recall some basic notions and 
results on nonclassical (conditional) symmetries of partial differential equations. 
This will form the basis for our discussion of myths in the next sections. 

Consider an involutive family Q = {Q 1 , . . . , Q } of I (I ^ n) first order differ- 
ential operators (vector fields) 

Q s = C(x, u)di + rj s (x, u)d u , s = l,...,l, 

in the space of the variables x and u, satisfying the condition rank ||£ Si (x, u)\\ = I. 

Here and in what follows x denotes the n-tuple of independent variables (x\, . . . , 
x n ), n > 1, and u is treated as the unknown function. The indices i and j run 
from 1 to n, the indices s and a run from 1 to I, and we use the summation con- 
vention for repeated indices. Subscripts of functions denote differentiation with 
respect to the corresponding variables, di = d/dxi and d u = d/du. Any function 
is considered as its zero-order derivative. All our considerations are in the local 
setting. 

The requirement of involution for the family Q means that the commutator of 
any pair of operators from Q belongs to the span of Q over the ring of smooth 
functions of the variables x and u, i.e., 

Vs,s' 3( ss,<7 = ( ss '°(x,u): [Q s ,Q s ']=C ss ' a Q' J . 

The set of such families will be denoted by £l l . 

Consider an rth-order differential equation C of the form L[u] := L(x, u/ r \) = 
for the unknown function u of the independent variables x. Here, U( r ) denotes 
the set of all derivatives of the function u with respect to x of order not greater 
than r, including u as the derivative of order zero. Within the local approach the 
equation C is treated as an algebraic equation in the jet space J r of the order r 
and is identified with the manifold of its solutions in J r . Denote this manifold 
by the same symbol C and the manifold defined by the set of all the differential 
consequences of the characteristic system Q[u] =0 in J r by Q( r ), i.e., 

Q (r) = {(x,n (r) ) G J r \D^ ...D^Q s [u] = 0, a { G NU{0}, [a| < r}, 

where Di = d Xi + u a+ $ t d Uci is the operator of total differentiation with respect 
to the variable Xi, Q s [u] := r] s — ^ m Ui is the characteristic of the operator Q, 
a = (ai, . . . , a n ) is an arbitrary multi-index, \a\ := a\ + • • • + a n , 5i is the 
multiindex whose ith entry equals 1 and whose other entries are zero. The variable 
u a of the jet space J r corresponds to the derivative d^u/dx" 1 . . . dx" n . 
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Definition 1. The differential equation C is called conditionally invariant with 
respect to the involutive family Q if the relation 

Ql r) L(x,u {r) )\ CnQ ^ = (1) 

holds, which is called the conditional invariance criterion. Then Q is called an in- 
volutive family of conditional symmetry (or Q-conditional symmetry, nonclassical 
symmetry, etc.) operators of the equation C. 

Here the symbol Qf r \ stands for the standard rth prolongation of the opera- 
tor Q s [IS E9]: 

Q(r)=Q S + E [D? 1 ...D^Q s [u]+e i U a +8 i )d Ua . 

0<|a|^r 

3 Myths on name and definition 

We restrict our consideration mainly to the case of families consisting of single 
operators (Z = 1) for simplicity and since mostly this case is investigated in the 
literature. Then the involution condition degenerates to an identity and we can 
omit the words "involutive family" and talk only about operators. 

Myth 1. A nonclassical symmetry operator Q of an equation C is a vector field Q 
which is a Lie symmetry operator of the united system of the equation C and the 
invariant surface condition Q[u] =0 corresponding to Q. 

This is the conventional non-rigorous way in order to quickly define nonclas- 
sical symmetry (see, e.g., |10t lllj). It becomes rigorous only after a special in- 
terpretation of the notions of system of differential equations and Lie symmetry. 
Otherwise, using the empiric definition leads to a number of inconsistencies. 

A closer look reveals that the above definition is a tautology. Indeed, the 
invariant surface condition Q[u] = means that the function u is a fixed point 
of the one-parametric local group Gq of local transformations generated by the 
operator Q. Therefore, we can reformulate the definition in the following way. 

Reformulation. If the set of those solutions of the equation C which are fixed 
points of Gq, is invariant with respect to Gq, then Q is called a nonclassical 
symmetry operator Q of the equation C. 

The tautology of the reformulation is obvious. If each element of the set is 
invariant then the whole set is necessarily invariant. The definition of nonclassi- 
cal symmetry according to Myth Q] leads to the conclusion that any differential 
equation is invariant, in the nonclassical sense, with respect to any vector field in 
the corresponding space of dependent and independent variables. 

The case when the equation C has no Q-invariant solutions fits well into the 
non-rigorous approach in the sense that the empty set is a particularly symmet- 
ric set. 
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Therefore, uncritically following the non-rigorous approach, we would get no ef- 
fective methods for constructing exact solutions and no information on the partial 
differential equations under consideration. 

There exist a number reformulations of Myth Q] in the literature in different 
terms. The first one is in terms of conditional symmetry. 

Myth 2. A nonclassical symmetry operator Q of an equation C is a conditional 
symmetry operator of the equation C under the auxiliary condition Q[u] =0. 

The association of nonclassical symmetries (under the name Q-conditional sym- 
metries) with conditional ones can be traced back to [7] (see also [8] and earlier 
papers of the same authors). Here the term conditional symmetry is understood 
in the following sense [6] (it can easily be defined for the case of a general system 
of differential equations). 

Definition 2. A vector field Q is called a conditional symmetry operator of a 
system C of differential equations under an auxiliary condition £ (which is another 
system of differential equations in the same variables) if Q is a Lie symmetry 
operator of the united system of C and £ . 

Conditional symmetries defined in this way essentially differ from nonclassi- 
cal symmetries. In particular, auxiliary conditions for conditional symmetries 
do not involve any associated conditional symmetry operators. The conditional 
symmetry operators of a system C under an auxiliary condition £' form a Lie 
algebra. Conditional symmetry indeed is a kind of symmetry and can be applied 
to generate new solutions from known ones. At the same time, in contrast to 
the case of nonclassical symmetries, finding auxiliary conditions associated with 
nontrivial conditional symmetries is an art rather than an algorithmic procedure. 
This is why sometimes nonclassical symmetries are called either Q -conditional 
symmetries, where the prefix "Q" is used to emphasize the differences between 
nonclassical and conditional symmetries, or conditional symmetries without any 
connection with Definition [2j 

The second reformulation of Myth [T] is in infinitesimal terms. Note that in- 
finitesimal criteria lie at the basis of Lie symmetry theory since they allow one 
to study linear problems for infinitesimal transformations instead of nonlinear 
problems for finite transformations. 

Myth 3. The conditional invariance criterion for an equation C and an oper- 
ator Q coincides with the infinitesimal Lie invariance criterion for the united 
system {C,Q[u] = 0} with respect to the same operator, i.e., 

Q( r )L[u] = if L[u] = and Q[u] = 0. 

The infinitesimal Lie invariance criterion for the invariant surface condition 
Q[u] =0 with respect to the operator Q is identically satisfied as an algebraic 
consequence of this condition since 



Q(r)Q[u] = Q{l)Q[u] = (Vu ~ tiuj)Q[u] = if Q[u] = 0. 
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We also have 

Q(r)L[u] = ?DiL[u] + Y, LuMD? 1 • ..D% n Q[u], (2) 

\a\^r 

i.e., the equation Q^L[u] = is a differential consequence of the equations 
L[u] = and Q[u] = and, therefore, becomes an identity on the set of their 
common solutions. This tautology was first observed in |17j . 

In the local approach to group analysis of differential equations, a system 
of differential equations is associated with the infinite tuple of systems of alge- 
braic equations defined by this system and its differential consequences in the 
infinite tower of the corresponding jet spaces. The exclusion of the differential 
consequence Qr r \L[u] when considering the system L[u] = and Q[u] = seems 
unnatural from the viewpoint of group analysis. 

A variation of Myth [3] is to replace, due to the Hadamard lemma, the "invari- 
ance condition" holding on the solution set of the system L[u] =0 and Q[u] = 
by the associated multiplier-condition, to be satisfied on the entire jet space J r . 

Myth 4. An operator Q is a nonclassical symmetry of an equation C if there 
exist A 1 and A 2 such that 

Q ( , r) L[u] = \ 1 L[u] + \ 2 Q[u}. (3) 

The problem is to precisely define the nature of the multipliers A 1 and A 2 . 
A number of different conditions on the multipliers have been put forward in the 
literature. The simplest version is to prescribe no conditions at all on A 1 and A 2 , 
which is obviously unacceptable. 

Sometimes A 1 and A 2 are assumed to be differential functions. This condition 
is natural for A 1 but overly restrictive for A 2 . In fact, if only such A 2 are allowed, 
the equivalence relation of nonclassical symmetries up to nonvanishing functional 
multipliers will be broken. Moreover, in this case the associated invariance cri- 
terion will become merely a sufficient condition for an ansatz constructed with 
the operator Q to reduce the equation C. As a result, a number of well-defined 
reductions may be lost. 

On the other hand, requiring that both the multipliers A 1 and A 2 are polyno- 
mials of total differentiation operators with respect to the independent variables, 
whose coefficients are differential functions, is too weak an assumption. It arises 
from the association of nonclassical symmetries with conditional symmetries for 
which such multipliers are admissible. If we choose 

A 1 = e Di and A 2 = ]T L Ua [u]B? . . . £>£» , 

\a\^r 

condition ([3]) obviously becomes an identity for any operator Q. In other words, 
condition ([3]) reduces to the tautology (|2|) if both A 1 and A 2 are treated as differ- 
ential operators of the above kind. 
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Comparing Definition Q] and Myth 0] shows that A 1 should be a differential 
function (i.e., a zeroth order operator) and A 2 should be an order (r — 1) operator. 
These conditions for the multipliers can be weakened. Thus, bounding the order 
of total differentiations in A 2 is not essential. If A 1 is a differential function, 
condition ([3]) implies that A 2 cannot include total differentiations of orders greater 
than r — 1. At the same time, explicitly prescribing the bound allows one to fix 
the order of the jet space under consideration. 

Myth 5 (The main myth of the theory). Nonclassical symmetry is a kind of 
symmetry of differential equations. 

Any kind of symmetry of differential equations (Lie, contact, hidden, con- 
ditional, approximate, generalized, potential, nonlocal etc.) has the invariance 
property, i.e., symmetries transform solutions to solutions in an appropriate sense. 

The basic prerequisite of the definition of nonclassical symmetry is the consid- 
eration of only the set of solutions invariant under the associated finite transfor- 
mations. It is impossible to use nonclassical symmetries in order to generate new 
solutions from known ones. A nonclassical symmetry operator Q of C represents 
only a symmetry of 

• each Q-invariant solution of C (as a weak symmetry |17| ) and 

• the manifold C n Qm in J r , where r = ord£. 

The manifold CnQ^ is properly related to the joint system L[u] = and Q[u] = 
of differential equations only if the operator Q and the equation C satisfy the 
conditional invariance criterion. 

At the same time, properties of the set of nonclassical symmetries and proper- 
ties of the set of Q-invariant solutions for each nonclassical symmetry operator Q 
characterize the equation C 

Since a nonclassical symmetry is not in fact a kind of symmetry of differential 
equations, it is of utmost importance to discuss possibilities for replacing the name 
by one not involving the word "symmetry". 

4 Nonclassical symmetry, compatibility and reduction 

To understand the real nature of nonclassical symmetry, we discuss properties 
and applications of Lie symmetries and single out those of them which carry over 
to nonclassical symmetries. 

Properties of Lie symmetries: 

Invariance. Any Lie symmetry (in the form of a parameterized family of finite 
transformations) locally maps the solution set of the corresponding system of 
differential equations onto itself. This is the main characteristic of any kind 
of symmetry. It gives rise to the possibility of generating new solutions from 
known ones. 
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Formal compatibility. Attaching the invariant surface conditions associated with 
a Lie invariance algebra to the initial system of differential equations results 
in a system having no nontrivial differential consequences. In other words, 
the invariant surface conditions forms a class of proper universal differential 
constraints and, therefore, is appropriate for finding subsets of solutions of the 
initial system. 

Reduction. Each Lie invariance algebra satisfying the infinitesimal transversality 
condition leads to an ansatz reducing the initial system to a system with a 
smaller number of independent variables, i.e., the reduced system is more easily 
solvable than the initial one. 

Conditional compatibility. There exists a bijection between solutions of the ini- 
tial system which satisfy the invariant surface conditions, and solutions of the 
corresponding reduced system. This means that all solutions of the initial sys- 
tem invariant with respect to a Lie invariance algebra, can be constructed via 
solving the corresponding reduced system. 

For nonclassical symmetries, the property of invariance is broken but the other 
properties (formal compatibility, reduction, conditional compatibility) are pre- 
served. In fact, the conditional invariance criterion (JTJ) is the condition of formal 
compatibility of the joint system L[u] = and Q[u] = [23]. We can identify 
nonclassical symmetries of C with first- order quasilinear differential constraints 
which are formally compatible with £0 

Definition 3. The differential equation C is called conditionally invariant with 
respect to the involutive family of operators Q if the joint system of C with the 
characteristic system Q[u] = is formally compatible. 

1 ln fact, this claim and Definition [3] are not entirely rigorous. Their precise formulation 
depends on what definition of formal compatibility is used. Consider, e.g., the definition presented 
in |24II25| . We temporarily use notations compatible with these references, hence slightly different 
from the rest of the paper. 

Let C r be a system of I differential equations L 1 [it] =0, . . . , L l [u] = in n independent vari- 
ables x = (;ri, . . . , x n ) and m dependent variables u = (it , . . . , u m ), which involves derivatives 
of u up to order r. The system £ r is interpreted as a system of algebraic equations in the 
jet space J r and defines a manifold in J r , which is also denoted by L r . The sth order pro- 
longation C r +s of the system £ r , s £ N, is the system in J r+S consisting of the equations 
D i 1 . . . D"" L [u] = 0, k = 1, . . . , I, \a\ ^ s. The projection of the corresponding manifold on 
J r+S ~ q , where q £ N and q ^ s, is denoted by C^_ g _ q . The system L r is called formally integrable 
(or formally compatible) if £^+s ~ C-r+s for any s £ N |24l 125] . 

The first obstacle in the harmonization of the above definition of formal compatibility and 
the definition of nonclassical symmetry is that the equations L[u] =0 and Q[u] — have, as 
a rule, different orders. Therefore, differential consequences of the equation Q[u] — should 
be attached to the joint system L[u] =0 and Q[u] =0 before testing its compatibility. The 
second obstacle is that the order of L[u] may be lowered on the manifold Qr r ) if Q is a singular 
vector field for the equation L[u] = 0. Hence instead of the equation L[u] =0 we should use 
the equation L*[u] = 0, where L» is a differential function which coincides with XL on Q( r ) 
for some nonvanishing differential function A and whose order r* is minimal among differential 
functions possessing this property. Finally, we arrive at the following definition: The differential 
equation C is called conditionally invariant with respect to the involutive family of operators Q 
if the system L t [u] = 0, D" 1 . . .D" n Q[u] = 0, \a\ < r», is formally compatible. 
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What is the main property that adequately represents the essence of nonclas- 
sical symmetry? 

The fact that the characteristic equations Q s [u] = are quasilinear and of first 
order implies the possibility of integrating them explicitly, i.e., an ansatz associ- 
ated with the characteristic system Q[u] =0 can be constructed. In view of the 
Frobenius theorem, the involution and transversality conditions for the family Q 
(together with the fact that the operators from Q are of first order) imply that 
the ansatz involves one new unknown function of n — I new independent vari- 
ables. Then the formal compatibility of the joint system L[u] = and Q[u] = 
guaranties the reduction of C by the ansatz to a single differential equations £ 
in n — I independent variables. Thus, the number of dependent variables and 
equations are preserved under the reduction with Q and the number of indepen- 
dent variables decreases by the cardinality of Q, i.e., similarly to Lie symmetries 
nonclassical symmetries lead to the conventional reduction of the number of in- 
dependent variables. 

There exist integrable differential constraints which are not formally compat- 
ible with the initial system. Differential constraints can be formally compatible 
with the initial system and, at the same time, non-integrable in an explicit form. 
An ansatz constructed with a general integrable differential constraint may involve 
a number of new unknown functions depending on different variables. Therefore, 
only all the above properties combined (first order, quasilinearity, formal compat- 
ibility, transversality and involution) result in the classical reduction procedure!! 

The conditional invariance of the equation £ with respect to the family Q is 
equivalent to the ansatz constructed with this family reducing C to a differential 
equation with n — I independent variables |28j . Moreover, reducing the number 
of independent variables in partial differential equations is the main goal in the 
study of nonclassical symmetries. Since the reduction by the associated ansatz 
is the quintessence of nonclassical symmetries, it was proposed in [211 E21 ES] to 
call involutive families of nonclassical symmetry operators families of reduction 
operators of C. 

Another important property holding for Lie symmetries is broken for nonclas- 
sical symmetries. Let the equation C be of order r and 

L (k) = {D? . . . D^L[u) = 0, \a\^k- r}. 

Denote by Lqa a maximal set of algebraically independent differential conse- 
quences of C that have, as differential equations, orders not greater than k. We 
identify L(u\ with the corresponding system of algebraic equations in J k (x\u) and 

2 Extended notions of reduction are also used. Thus, weak symmetries imply reductions de- 
creasing the number of independent variables, preserving the number of unknown functions and 
increasing the number of equations [17]. The reduced system can be much more overdeter- 
mined than the initial one. The reductions associated with higher-order nonclassical symmetries 
preserve the determinacy type of systems, simultaneously increasing the numbers of unknown 
functions and equations |15j . 
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associate it with the manifold Cn<.\ determined by this system. For Lie symmetries 
we have the following properties. 

1. If Q is a Lie symmetry operator of Cf r \ then Q is a Lie symmetry operator of 

for any p > r. 

2. If Q is a Lie symmetry operator of £( p ) for some p > r then Q is a Lie symmetry 
of £ (r) . 

The first of these properties extends to nonclassical symmetries but this is not 
the case for the second one. In fact: 

1. If Q is a Lie symmetry operator of Ct r \ n Qr r ) then Q is a Lie symmetry operator 
of £( p ) n Q (p ) for any p > r. 

2. The fact that Q is a Lie symmetry operator of £( p ) n Q( p \ for some p > r does 
not imply that Cr r ) H Q( r ) admits the operator Q. 

Example 1. Let L[u] = + u xx + tu^, £: L[u] = and Q = dt- Then the 
manifold £( 2 ) ^ 2(2) is determined in J 2 by the equations ut = Utt = ut x = and 
u X x = —tu x . Since 

the operator dt is not a reduction operator of C. Substituting the corresponding 
ansatz u = <p(uj), where the invariant independent variable is u = x, into C results 
in the equation {p WUJ + tip u = 0, in which the "parametric" variable t cannot be 
excluded via multiplying by a nonvanishing differential function. As expected, 
the ansatz does not reduce the equation C 

Consider the same operator Q and the first prolongation £/ 3 \ of C, which is 
determined by the equations L[u] =0, D t L[u] = and D x L[u] = 0. The manifold 
£(3) H Q(3) is singled out from J 3 by the equations 

U t = U tt = U tx = Uut = Uttx = Utxx = 0, U x = U xx = U XX x = 0. 

The conditional invariance criterion is satisfied for the prolonged system Cr%\ and 
the operator Q: 

^(2) L l£ (3) ns (3) = Q W DtL \c (s) n Q{3) = Q(V D * L \c (s) nQ {3) = °> 

i.e., Q is a nonclassical symmetry operator of the system Crg\ and the above 
ansatz reduces £( 3 ) to the system of three ordinary differential equations (fu, = 0, 
^ = and (puLou = since 





and 



t 1 
1 
t 1 



^0. 



Note 1. In general, for any system C and any involutive family Q there exists 
an order r such that Ct r ) H Qm is invariant with respect to Qi r y This gives the 
theoretical background of the notion of weak symmetry |17j . 
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5 Definition of nonclassical symmetries for systems 

Myth 6. The definition of nonclassical symmetry for systems of differential equa- 
tions is a simple extension of the definition of nonclassical symmetry for single 
partial differential equations to the case of systems. 

Example [JJ and Note Q] indicate problems arising in attempts of defining non- 
classical symmetries for systems of partial differential equations. 

Let C denote a system L(x,um) = of / differential equations L 1 = 0, . . . , 
L — for m unknown functions u = (u 1 , . . . ,u m ) of n independent variables 
x = (xi, . . . , x n ). It is always assumed that the set of differential equations form- 
ing the system under consideration canonically represents this system and is min- 
imal. The minimality of a set of equations means that no equation from this set 
is a differential consequence of the other equations. By Lt^\ we will denote a 
maximal set of algebraically independent differential consequences of C that have, 
as differential equations, orders not greater than k. We identify Lr^) with the 
corresponding system of algebraic equations in the jet space J k and associate it 
with the manifold Ci k \ determined by this system. Let Lm = {L u , u = 1, . . . , I}. 
Note that the general system includes equations of different orders. 

What is the correct conditional invariance criterion for the system £? 

Q(r) L "\c {r) n Q(r) =^ M = 1 >---> Z? 
^M^1/: (r) ns (r) =0 ' " = 1 >-..J? 

All of the above candidates for the criterion are not satisfactory. The second 
candidate is not a good choice since it neglects the equations having lower or- 
ders than the order of the whole system. Taking the third candidate, we obtain 
nonclassical symmetries of a prolongation of the system. As shown by Exam- 
ple [H these may be weakly related to nonclassical symmetries of the system. It 
is not well understood what differential consequences are really essential. Thus, 
elements of Cm whose trivial differential consequences also belong to Cm are 
neglected by this candidate. 

Although all operators satisfying the first of the above criteria give proper 
reductions, it is overly restrictive and in fact is only a sufficient condition for 
nonclassical symmetries. Even Lie symmetries can be lost when employing it. 

The above discussion is illustrated by the following example. 

Example 2. Consider the system 

ut + (u ■ V)u — An + Vp + x x Vdivu = 0, divn = 0. (4) 

which is obviously equivalent to the system of Navier-Stokes equations describing 
the motion of an incompressible fluid. (The additional term x x V(divu) van- 
ishes if div-u = 0.) If we do not take into account differential consequences of 



12 



Michael Kunzinger and Roman O. Popovych 



system we derive the unnatural claim that this system is not conditionally 
invariant with respect to translations of the space variables X{. At the same time, 
the infinitesimal generators of these translations belong to the maximal Lie invari- 
ance algebra of the Navier-Stokes equations. A maximal set Lm\ of algebraically 
independent differential consequences of C that have, as differential equations, 
orders not greater than 2 is formed by the equations 

u t + {u- V)u — Au + Vp = 0, div u = 0, 
div u t = 0, V div u = 0, u^uj + Ap = 0. 

Here the indices i and j run from 1 to 3. The equation Q(2) div u = is identically 
satisfied on the set £( 2 ) n Q(2)- Therefore, the application of the second or third 
candidate for the conditional invariance criterion to the equation div u = gives 
no determining equations for nonclassical symmetries of the system 

Definition [3] can also not be extended to the case of systems in an easy way. 
The problem again is to define what set of differential consequences of the initial 
system should be chosen for testing formal compatibility with the appropriate 
characteristic system. 

The notions of nonclassical symmetry and reduction are strongly related in the 
case of single partial differential equations. It therefore seems natural for these 
notions to also be closely related in the case of systems. Hence the problem of 
rigorously defining nonclassical symmetries for systems is additionally complicated 
by the absence of a canonical extension of the classical reduction to the case 
of systems. A chain of simple examples can be presented to illustrate possible 
features of such an extension. 



6 Myths on number of nonclassical symmetries 

Myth 7. The number of nonclassical symmetries is essentially greater than the 
number of classical symmetries. 

At first sight this statement seems obviously true. There exist classes of par- 
tial differential equations whose maximal Lie invariance algebra is zero and which 
admit large sets of reduction operators. This is the case, e.g., for general (1 + 1)- 
dimensional evolution equations. At the same time, certain circumstances signif- 
icantly reduce the number of essential nonclassical symmetries. We briefly list 
them below. 

• The usual equivalence of families of reduction operators. Involutive fami- 
lies Q and Q of I operators are called equivalent if Q s = \ SCT Q a for some 
X sa _ x«r( x ,u) with det ||A SCT || + 0. 

• Nonclassical symmetries equivalent to Lie symmetries. 
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• The equivalence of nonclassical symmetries with respect to Lie symmetry 
groups of single differential equations [131 I2Q] and equivalence groups of 
classes of such equations [22] . 

• No-go cases. The problem of finding certain wide subsets of reduction oper- 
ators may turn out to be equivalent to solving the initial equation [12\ [2Tj . 

• Non-Lie reductions leading to Lie invariant solutions. 

Thus, the existence of a wide Lie symmetry group for a partial differential 
equation C complicates, in a certain sense, finding nonclassical symmetries of C 
Indeed, any subalgebra of the corresponding maximal Lie invariance algebra, sat- 
isfying the transversality condition, generates a class of equivalent Lie families of 
reduction operators. If a non-Lie family of reduction operators exists, the action of 
symmetry transformations on it results in a series of non-Lie families of reduction 
operators, which are inequivalent in the usual sense. Therefore, for any fixed value 
of I the system of determining equations for the coefficients of operators from the 
set O 1 (£) of families of I reduction operators is not sufficiently overdetermined to 
be completely integrated in an easy way, even after factorizing with respect to the 
equivalence relation in £r(£). To produce essentially different non-Lie reductions, 
one has to exclude the solutions of the determining equations which give Lie fam- 
ilies of reduction operators and non-Lie families which are equivalent to others 
with respect to the Lie symmetry group of C. As a result, the ratio of efficiency 
of such reductions to the expended efforts can become vanishingly small. 

7 Conclusion 

Although the name "nonclassical symmetry" and other analogous names for re- 
duction operators, which refer to symmetry or invariance, do not reflect actual 
properties of these objects, the usage of such names is justified by historical con- 
ventions and additionally supported by the terminology of related fields of group 
analysis of differential equations. It is a quite common situation for different fields 
of human activity that a modifier completely changes the meaning of the initial 
notion (think of terms like "negative growth", "military intelligence", etc.). Em- 
piric definitions of nonclassical symmetry can be used in a consistent way if all 
involved terms and notions are properly interpreted. Nevertheless, as we have 
argued, the term reduction operator more adequately captures the underlying 
mathematical content. 

In this paper we discussed certain basic myths of the theory of nonclassical 
symmetries, pertaining to different versions of their definition and the estimation 
of their cardinality. Over and above these, there are a number of more sophis- 
ticated myths concerning, among others, the factorization of sets of nonclassical 
symmetries, involutive families of reduction operators in the multidimensional 
case, and singular sets of reduction operators. A discussion of such myths re- 
quires a careful theoretical analysis substantiated by nontrivial examples and will 
be the subject of a forthcoming paper. 
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